Defect mediated melting and the breaking of quantum double symmetries 
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In this paper, we apply the method of breaking quantum double symmetries" 1 " 1 ^ to some cases of 
defect mediated melting. The formalism allows for a systematic classification of possible defect 
condensates and the subsequent confinement and/or liberation of other degrees of freedom. We also 
show that the breaking of a double symmetry may well involve a (partial) restoration of an original 
symmetry. A detailed analysis of a number of simple but representative examples is given, where 
we focus on systems with global internal and external (space) symmetries. We start by rephrasing 
some of the well known cases involving an Abelian defect condensate, such as the Kosterlitz-Thouless 
transition and one- dimensional melting, in our language. Then we proceed to the non-Abelian case 
of a hexagonal crystal, where the hexatic phasA— is realized if translational defects condense in a 
particular rotationally invariant state. Other conceivable phases are also described in our framework. 

PACS numbers: 02.20.Uw 64.60.-i 61.30.-v 61.72.-y 



I. INTRODUCTION 

In two-dimensional (quantum) physics certain quan- 
tum groups and Hopf algebras play an important role 
because these extended symmetries allow for a treatment 
of ordinary and topological quantum numbers on equal 
footing, whereas the natural appearance of an R-matrix 
accounts for the topological interactions in the system^. 
This means that both topological defects and ordinary 
excitations (such as Goldstone modes) may appear in 
the same representation of the underlying - often hid- 
den - Hopf algebra A, even in non-Abelian situations 
where the mutual dependence of these dual quantum 
numbers would otherwise be untractable. For a brief, 
tailor made introduction to Hopf algebras and the no- 
tational conventions we adopt in this paper, we refer to 
Appendix A of a related paper—, or to a more general 
review!. These extended symmetry concepts have found 
interesting applications in the domain of Quantum Hall 
liquids-, crystals and liquid crystals^, (2+l)-dimensional 
gravityiS, and potentially also to rotating Bose-Einstein 
condensatesii*i£ and other systems described by confor- 
mal field theory . 

Once this type of (hidden) symmetry was identified, 
a formalism for the breaking of quantum double sym- 
metries was proposed^, by assuming the formation of 
condensates of ordinary (electric), defect (magnetic), or 
mixed (dyonic) type. As was expected, the ordinary con- 
densates reproduce the conventional theory of symmetry 
breaking, though the systematic analysis of the subse- 
quent confinement of topological degrees of freedom us- 
ing the braid group is not standard. In the papers just 
mentioned it was shown that the usual picture of symme- 
try breaking had to be augmented with significant novel 



ingredients. In the case of Hopf double symmetry break- 
ing we assume the condensate to be represented by a 
fixed vector \<pa) in some nontrivial representation of 
the Hopf algebra A. This leads to the definition of an 
intermediate algebra T as the suitably defined stabilizer 
subalgebra of the vacuum state: T was originally defined 
as the largest Hopf subalgebra of A whose elements a 
satisfy 



n«(a) | <M = e(o) | <t> ) V aeT 



(1) 
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where e(a) denotes the counit or trivial representation 
of the Hopf algebra in question. This analysis is fur- 
ther generalized in a related paper—, where it is argued 
that sometimes more general (non-Hopf) algebras play 
an essential role. We then have to distinguish between 
a right or a left intermediate algebra, T r and 7/ respec- 
tively, which are all related in the sense that T C T T C\%. 
The complication that arises at this " intermediate" level 
is that certain representations of T may have nontrivial 
braiding properties with the vacuum state. This in turn 
means that the vacuum cannot be single valued around 
a particle carrying that representation. If this happens 
to be the case, it has the physical implication that such 
particles (representations) have to be confined. The low 
energy theory of non-confined degrees of freedom is then 
characterized by yet a different algebra called U. So the 
breaking of Hopf symmetries involves three algebras: the 
unbroken algebra A, the intermediate algebra (be it T, T r 
or %), and the unconfined algebra U. The generic break- 
ing scheme is depicted in Figure ^ With this formalism 
of Hopf symmetry breaking available it is important to 
take a closer look at suitable physical situations to check 
whether the theory produces correct and useful results. 
This is the subject of the remainder of this paper. Some 
applications to non-Abelian nematic crystals are the sub- 
ject of another paper appearing in parallel 9 . 
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FIG. 1: A schematic of the quantum symmetry breaking for- 
malism. One has to distinguish three levels: The unbroken 
Hopf algebra D(H), the intermediate algebra T r and the effec- 
tive unconfined algebra U. This scheme should be compared 
with the simple scheme of ordinary symmetry breaking where 
one only has to distinguish between the unbroken group G 
and some residual symmetry group H C G. 



II. ABELIAN DEFECT CONDENSATES 

In this section we will treat some examples of 
"Abelian" condensates. The rather well known exam- 
ples of Abelian defect condensates are interesting because 
they cause confinement but at the same time an analy- 
sis of the unconfined particles reveals that a so called 
symmetry restoration takes place. In other words Hopf 
symmetry breaking by defects may lead to the restoration 
of ordinary symmetries. 



A. Kosterlitz-Thouless revisited 

The most well known example of a defect condensate is 
the one that mediates the Kosterlitz-Thouless transition 
in the 2-dimensional XY models. In the isotropic phase, 
the internal symmetry group is 50(2). Representations 
of SO{2) are labelled by an integer: n £ Z. For reasons 
that will become clear shortly, we also introduce parti- 
cles that transform under projective unitary representa- 
tions of S0(2), which is equivalent to saying that they 
transform under faithful representations of the universal 
covering group R of SO{2). These R representations are 
labelled by a real number A £ R, giving the irrep py. 



Px{x) 



x e 



(2) 



Thinking of A as labelling a spin degree of freedom one 
has that the particle is bosonic if A 6 Z, because p\(2tt) = 
1. If A eZ+|, the irrep is fermionic, because p\{2 j n) = 
— 1. For any other value of A, the irrep would be called 
an Abelian anyon. 

In the low temperature phase of our system, the SO(2) 
symmetry is completely broken by a bosonic (A = 1) 
condensate implying that the covering group R is broken 



to Z. The residual symmetry transformations, labelled 
by n £ Z, correspond to rotations by 27rn. The Hopf 
symmetry in this phase is the quantum double of Z de- 
noted by -D(Z), because the defects are labelled by an 
integer n, as a consequence of the fact that Hi(R/Z) = 
no(Z) = Z. The electric irreps in this phase are labelled 
by p £ [0, 27r] , corresponding to the irreps of Z. However, 
we can still include all the electric irreps of the isotropic 
theory. In other words, we can restrict the irreps of R to 
Z, giving representations of Z, some of which are equiva- 
lent. We like to take all the electric irreps along, because 
we want to see which irreps of the isotropic phase sur- 
vive unconfined in the various phases. Note that the 
Goldstone boson transforms trivially under Z, because 
27T rotations don't affect it. So in this phase we have the 
representations n™ of D (Z), where 



£>(Z) = F(Z) x CZ. 



(3) 



If we heat the system in this -D(Z) phase, the Kosterlitz- 
Thouless transition takes place corresponding to the for- 
mation of a condensate of the defect \4>q) = \n = 1). To 
see what happens one first has to determine the inter- 
mediate algebra T and then use the braid relations to 
determine which representations are confined to obtain 
the effective unconfined algebra U. The present case cor- 
responds to a single defect condensate and the general 
structure of T and U have been given elsewhere^. The 
outcome is 



% = F(H/(g A )) x CN A 

U = F(N A /(g A )) x CN A /(g A ). 



(4) 
(5) 



where (g A ) is the group generated by the group element 
g A corresponding to the defect in the condensate, and 
this group (g A ) is a normal subgroup of N A , which is the 
normalizer in H of the element g A . In the case we are 
considering, (g A ) — N A , so for U we obtain the rather 
trivial result: 

U = D(Z/Z) = D(e), 

from which we conclude that there are no surviving de- 
fects in the unconfined phase and also that from the set 
electric reps of H = Z (i.e. p £ [0, 2tt]) non are left. 

But what about the original set of spin representations 
Pjl? Considering the braiding of a representation with 
spin A with a defect \n), we obtain thatSi 



\n)\\)=j^\n)\\). 



(6) 



The only irreps of R that braid trivially with the conden- 
sate |1) are those labelled by an integer m £ Z: 



p m {x) 



x £ 



So the residual theory has no defects, and the unconfined 
electric irreps are labelled by n £ Z, in other words the 
fractional charges are confined. A theory with no defects 
and charges labelled by n £ Z is an 5*0(2) theory. This 
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argument suggests that the symmetry has been restored 
to a full 50(2), i.e. the original R theory has been com- 
pactified, and the "spin" is now quantized. 

In general, condensing the \n) defect gives a theory 
with Z„ defects (defects defined modulo n), and the elec- 
tric irreps are labelled by — , k £ Z, meaning that the 
symmetry is restored to SO(2)fZ n . 



B. 1-Dimensional melting 

In the previous example we studied a phase with in- 
ternal symmetry, we now turn to a similar example 
with a spatial translational symmetry. Consider a one- 
dimensional system, say defined along the x-axis. In the 
high energy phase, the external (space) symmetry group 
of the system is the translation group R. As in the previ- 
ous example, the irreducible unitary representations (ir- 
reps) p\ of the external symmetry group R are again 
labelled by A <G R (see ©). We now break this R sym- 
metry to Z, by the formation of a one dimensional crys- 
tal, i.e. a 1-dimensional lattice. The lattice sites are at 
positions x = . . . , —2a, —a, 0, a, 2a, . . . We can label the 
lattice sites by elements of Z, where m £ Z corresponds 
to the lattice site at position x = ma. 

The irreps fl^ of the residual external symmetry group 
Z are labelled by p, £ [0,27r/a): 

Q^m) = e ^ ma , m £ Z. (7) 

p is the momentum of the irrep Q^. In the broken phase 
with broken continuous translation invariance phonons 
appear as Goldstone modes. These phonons have mo- 
menta k £ R. Furthermore the broken phase is also 
characterized by possible defects (dislocations) labelled 
by an integer n because IIi(R/Z) = IIo(Z) = Z. 

Now imagine that we condense the defect n). If we 
now consider the braiding of a phonon with momentum 
k with this defect, we obtain again that 

Ti 2 -\n)\k) =e lkna \n)\k). (8) 

This braiding is trivial if and only if e lkna — 1, which 
is the case if k is an integer multiple of — . Since 
k £ £?, we have n phonons that braid trivially with 
the condensate, namely the phonons |fc) with momentum 
k = 0, 2zl,2^l . ( n _i)2zr The defects \m),m £ Z are 
unconfined, but they are now defined modulo since 
\n) is condensed. 

Thus the condensation of \n) has led to a theory with 
Z n defects, and Z„ phonons and we end up with an un- 
confined algebra U = D(Z n ). 

In particular, if we take n = 1, i.e. we condense the 
|1) defect, we obtain a phase with trivial symmetry, in 
which there are no defects, and all phonons are confined. 

Let us now discuss melting. If we condense the defect 
|1), then the unconfined modes are those with momen- 
tum k = 0, — 2 — . 3— , . . .. So all phonons are confined 



(since a phonon cannot have exactly zero crystal momen- 
tum), and the unconfined translational irreps have mo- 
mentum k £ {— n, n £ Z}. Thus we have a theory with 
momentum modes labelled by n £ Z, and no defects. 
This is obviously a theory with a periodic translational 
U(l) symmetry. This is exactly the same argument as the 
one for the Kosterlitz-Thouless phase, except that there 
the Goldstone modes were not confined, yet the symme- 
try was restored to SO{2) . In this theory one could of 
course also consider situations in which various mutually 
compatible condensates (that trivially braid with each 
other) occur simultaneously, leading to yet other resid- 
ual symmetries and excitations. 

C. A defect condensate in a 2-dimensional smectic 

We may extend the previous reasoning to a particular 
2-dimensional system, where we consider only transla- 
tional degrees of freedom. In the high energy phase, the 
(translational) symmetry of the system is R 2 . The irreps 
P£ of R 2 are labelled by k = (k xi k y ) £ R 2 : 

p e (x,y) = S k '* +k **\ (9) 

We break this symmetry to R x Z, by condensing the 
mode P(o,a)- The system is then composed of a set of 
lines parallel to the x-axis. In the y-direction, we have 
the same situation as above, i.e. the symmetry is bro- 
ken from R to Z in the y-direction. In the x-direction, 
the symmetry is still unbroken, so we can neglect the 
x-direction. Let us be more precise about this. The 
phonons, as Goldstone modes of the broken translations 
in the y-direction have a momentum that points in the 
y-direction so they have k x = 0. 

Now condense the basic dislocation in the y-direction, 
the |1) defect. The phonons all braid nontrivially with 
the condensate! Just as above, all the phonons are con- 
fined. 

Now that we've seen how the symmetry is restored, we 
can also discuss what happens in the Abrikosov lattice in 
Type II superconductors, for example, in our terminol- 
ogy. The symmetry is first fully restored by the defect 
condensate, and then it is broken to a lattice, as it were 
by the condensation of two defect density modes, similar 
to the formation of the smectic. It is a dual version of 
what we just discussed. 

III. NON-ABELIAN DEFECT CONDENSATES 

We now want to move on to defect condensates in a 
non-Abelian setting, and it turns out that we then have 
to distinguish different types of defect condensates. Con- 
sider a phase described by a quantum double D(H), or a 
generalized quantum double F(H m ) ® CiJ e ;2i, and pick 
a purely magnetic representation IT^ (1 is the trivial rep- 
resentation of the centralizer Na)- A basis of the vector 
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space on which this irrep acts is given by {Iff/ 1 )}, where 
the gf are the different defects in the conjugacy class 22 
Ca whose representative element we denote by qa- In 
that defect class Ca we may consider the following in- 
equivalent types of condensates: 

• Single defect condensate 

= \gf) (io) 

• Class sum defect condensate 

<t> = E \at)=-\c A ). (ii) 

We denote the condensate by \Ca), where gA is the 
preferred element of 

• Combined defect condensate 

0=Z>> (12) 

where E is some subset of the defects in one class. 
We need only take the elements to be within one 
class because it turns out that we need only study 
the cases where the condensate is the sum of vectors 
in the same irrep. 

The single defect and class sum defect condensates are a 
special case of the combined defect condensate. In the 
following we will encounter both class sum and combined 
defect condensates. 



A. The hexagonal crystal and its regular modes 

We start with the planar, achiral, hexagonal crystal 
which has a space symmetry group Z 2 x Zq. Strictly 
speaking in this case the defect classes (i.e. the sets of 
defects transformed into each other under residual global 
symmetry transformations) are larger than the conjugacy 
classes of the magnetic group, but for our present pur- 
poses it suffices to restrict our considerations to the con- 
jugacy classes. 

The first Brillouin zone (rescaled) is shown in Fig0 it 
is a hexagon. Points in the Brillouin zone correspond to 
irreps of the translation group. Namely, a point in the 
Brillouin zone is a momentum vector k — (k x ,k y ), and 
the associated irrep of Z 2 is: 

(n, m ) h-> e*"^™'. (13) 

B. Topological defects in the hexagonal phase 

The homotopy groups which define the defects are: 

Il 1 (G/H) = H (14) 
TL 2 (G/H) = 0, (15) 



ir 



FIG. 2: A triangular lattice, representing a crystal with 
hexagonal symmetry. The dashed vectors are vectors in the 
reciprocal lattice, and the shaded hexagon represents the first 
Brillouin zone. 



where H is the lift of H in the universal covering group 
G of G. In three dimensions, there are no monopolcs 
because the space group H is discrete. The line defects 
are characterized by elements of H. 
First we consider a two-dimensional hexagonal crystal. 
To determine H we need to lift H to the covering group 
of G = K 2 x SO(2). The translational part M 2 is already 
simply connected. The covering group of SO(2) is K. An 
element 8 of M corresponds to a rotation over an angle 6. 
In the covering group, however, a rotation over 2n is no 
longer equivalent to the identity. Similarly, the covering 
of Z 6 is 

%, = Z x Z 6 ~ Z (16) 

Therefore Z^xZ 6 ~ Z 2 x Z, and 

Ui(G/H) ~ Z 2 x Z. (17) 

Therefore a defect is an element of Z 2 x Z. We denote it 
by (n, m)r p , with n,m,p G Z. r is a y rotation, (1,0) a 
translation by one lattice vector in the x-direction, and 
(0,1) a translation by one lattice vector in the y-direction. 
A defect of the form (n, m) (i.e. with p = 0) is a dislo- 
cation. An r p defect (n — m = 0) is a disclination. 
Note that the structure of defects described by (|17fl is in- 
deed non-Abelian, which means that we have to think of 
them forming multiplets corresponding to the conjugacy 
classes of the first homotopy group. We return to this 
question in the following section when we construct the 
Hopf algebra describing this phase. 

Note that although both rotational and translational 
Goldstone modes are present in a two-dimensional crys- 
tal, it is known that the rotational Goldstone mode, 
called the roton, is much more massive^* than the trans- 
lational modes, because it costs a lot of energy to rotate 
different parts of the crystal relative to each other. For 
the same reason, disclinations are much more massive 
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than dislocations, and at low temperatures we only have 
dislocations and translational Goldstone modes. 

For a three-dimensional crystal, the analysis is very 
similar. For example, in a crystal made up of layers of 
hexagonal crystal, G = K 3 x SO(3) is broken to H = Z 3 X 
Zq, so the defects have an extra label, which represents 
translational defects in the direction perpendicular to the 
layers: 



(18) 



For a dislocation characterized by (n, m,fc), the vector 
(n, m, k) is the Burgers vector of the defect. A general 
dislocation has both screw and edge character, and the 
Burgers vector is given by the integral of the order pa- 
rameter field u[x) along a loop L encircling the corei^. 
The fusion rules of line defects have some interesting fea- 
tures, which are crucial for the story to come. We illus- 
trate these features for the two-dimensional hexagonal 
crystal. To determine the fusion rules, we need the con- 
jugacy classes of the defects. This requires knowledge of 
the multiplication in Z 2 x Zg, which is set by 



r{a,b)r~ 1 = (-6, a + b). 



(19) 



The conjugacy classes are listed in Table [I] From this we 
can for example prove that two disclinations can fuse to 
form any dislocation in this crystal. Namely, the fusion 
rules are given by multiplying entire conjugacy classes. 
If we fuse C r and C r -i, we can get any dislocation 
(m,n) because {m,n)r is an element of C r . Thus, from 
a topological point of view a dislocation can "dissociate" 
into two disclinations. 



C. The double symmetry of the hexagonal crystal 

In the previous subsection we saw that the magnetic 
group of a hexagonal crystal with large symmetry trans- 
formations is Z 2 x Z. Let us now consider a hexago- 
nal crystal and include the reflections, then the magnetic 
group is the same^i: H m = Z 2 x Z. The electric group 

is H el = Z 2 ^Z) 6 = (Z 2 x D 6 ) x Z. The last Z part 
corresponds to rotations of 2im, n € Z, which we denote 



by 



,,6?! 



These rotations are in the center of the group. 



Thus the Hopf symmetry is 25 : 

A = F(Z 2 x Z) <g> C((Z 2 x D 6 ) x Z). 



(20) 



We write an element of (Z 2 x D 6 ) x Z as (k, l)s J r k+6n , 
with k,l,n G Z, j G Z2, n G Z. The reflection in Dq is 
denoted by s. We take s to be the reflection with respect 
to the x-axis. 

To analyze the quantum numbers of the excitation sec- 
tors of our theory, we proceed as usual and have to con- 
struct the representations of the algebra given abov o 1 ^ 3 . 
The representations of the generalized quantum doubles 
have a generic structure. There is an orbit or defect class 



Representative gA 


Conjugacy class Ca 


r 6k , k G Z 


C r e k = {r 6k } 


(a, bjr 
k G Z, a, b > 


/ ' 

Lj (a,b)r 6k — 

{(a,by k ,(~-b,a + b)r ak , 
(-a-b, a y' k ,(-a y -b)r 6k , 
(Oj —a — b)r , (a + o, — a)r j 


r 1+6fc , k G Z 


Cyl+6fc = 

{(m,n)r f : (m, n) £ Z } 


r - 1+6k , fc G Z 


C r -l+6fc = 

r — l+6fc . /™ „\ r- n?2\ 


r 2+6fe , k G Z 


C r 2 + 6fc = 

\ l?Tlj 71 )r . ?7l — 71 t o-^j 


(l,0)r 2+6fc , k G Z 


^(l,0)r 2 + S fc — 

|(m,njr : m — n f, 6L\ 


r 3+ak , k G Z 


C r 3 + 6k = 


r 3+ak , k G Z 


C(l,0)r 3 + 6ft — 
\{m,n)r : (m,nj £ z£i x z/fj 


r 4+6fc , k G Z 


C r 4+6fc 

{(m,n)r 4+Bfc : m - n G 3Z} 


(l,0)r 4+6fc , k G Z 


C(l,0)r 4 + 6fc = 

{(m,n)r 4+6fc : m - n £ 3Z} 



TABLE I: The classification of defects for the 2-dimensional 
hexagonal crystal in classes of Z 2 x Ze = Z 2 x Z. 



of the group in the function part labelling the defects, 
paired with a representation of the centralizer Na for the 
ordinary/ Goldstone modes. These representations will 
be constructed in this section and the results are listed 
in Table Hm We need the orbits O of ((Z 2 x Z 6 ) x Z) 
under the action of (Z 2 x D§) x Z, the centralizer N of a 
chosen preferred element of each orbit, and the irreps of 
the centralizer. The action of (Z 2 x D 6 ) x Z on Z 2 x Z 6 
is given by: 



r(a,6)r" 1 = (-b,a + b), 
s(a, 6)s _1 = (a + b, —b). 



srs 1 = r 1 
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Orbits 





& r 6k — C r 6k 




{C(a,b)r® k i C( a +6,-6)r 6fe } 


O r l + 6k 


= {C r l + 6fc , C r -l-6k } 


O r 2 + 6k 


— {C r 2 + 6fc , C r —2-Gk } 



C(l,0)r 2 + 6 fc — {C(l,0)r 2 + 6fc iC(l,0)r-2-8fc} 

C r 3 + 6fc = {C r 3 + 6/s , C r -3-6fc } 
®(l,0)r 3 + ek ~ {C(l,0)r3 + 6fc iC(l,0)r-3-6fc} 



TABLE II: The orbits of Z 2 x Z under the action of (Z 2 x 
Z)e) x Z. The conjugacy classes are defined in Table □ 



There is one class of orbits, O r 6k , for which determin- 
ing the irreps of the centralizer (Z 2 xi Dq) x Z is a little 
complicated. We can determine the irreps of Z 2 XI Dq 
and those of Z separately, because the irreps of a direct 
product of groups are the tensor products of the irreps 
of the separate groups. The irreps of Z 2 xi Dq can be 
constructed by the method of induced representations 
described in the appendix. We must pick a momentum 
vector k which corresponds to an irrep of the transla- 
tion group Z 2 , and act on it with Dq . The elements that 
transform k into a vector that corresponds to the same 
irrep form a subgroup of Dq, which we denote by {Dq)^. 
Thus we are considering orbits in the first Brillouin zone 
(note that these classes have nothing to do with the or- 
bits listed in Table HJ). We must find all possible cases of 
vectors in the first Brillouin zone, with different (Z)g)^. 
For the hexagonal crystal, there are seven different types 
of orbits, and a representative element of each type is 
given in Figure |21 



Conjugacy classes 


N 


Irreps FI^ 


C r 6k 


(Z 2 -a As) x Z 


TT'' 6 " 


fcez 






C(a,b)r f >'= > C(a + b,-i))r6'= 


(Z 2 ) x Z 


u (a,b)r 6k 


a, b > 0; b > a; k G Z 




kern 2 




(Z 6 ) x 1 


1+6, 


k € Z 




m e Ze 


C r 2 + 6k , C r — 2 — 6k 


(Z 6 ) x Z 




k € Z 




m G Ze 


C(l,0)r 2 + 6 fc > C(l,0)r- 2 - 6fc 


(Z 3 ) x Z 


n (l,0)r 2 + 6fc 


k € Z 




n e z 3 


C r 3+Gk , C r -3-6fc 


(Z 6 ) x Z 


3 + 6fc 

IT 


ken 




m G Z6 




(Z 2 ) x Z 


n (l,0)r 3 + 6fc 


ken 




i e z 2 



TABLE III: The irreps of F((Z 2 x Z) <g> C(Z 2 x As) x Z. The 
index A G [0, 2-7r) corresponds to the irreps of the Z part of 
the centralizers, which represents the rotations by an angle 
2-nn. The index p in the first representation listed, refers to 
the irreps of Z 2 xi D 6 defined in Table HVI 



(-1,1) 



(-1,0) 




(0,-1 



FIG. 3: The first Brillouin zone of the planar hexagonal crys- 
tal. The gray area is a set of representatives of orbits under 
rotations and reflections. The numbered squares are repre- 
sentatives of the different types of orbits. The x- and y-axes 
have been rescaled. The unit distance on these rescaled axes 
corresponds to a distance of — , where a is the lattice spacing 
in the crystal. 



We then call Z 2 xi (Dq)^ the little group of k. We must 
determine the irreps of the little group, and then induce 
irreps of the whole group. This procedure gives all the 
irreps. which are given in Table IrTI 
As an example, We explicitly give the irrep corresponding 
to orbit 5 in Figure |3] In the appendix it is illustrated 
how to derive this. The matrices of the irrep are 



2 a^i 2 b 



(a,b) 











f o 




S I ► 









1 g^n 
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Orbit number 


(AOs 


Number of irreps 


Dim. of irreps 


1 




6 


1, 1, 1, 1, 2, 2 


2 


Z 2 


2 


6, 6 


3 


Z 2 


2 


6, 6 


4 


e 


1 


12 


5 


Z 2 x] Z 2 


4 


3 . 3 j 3 j 3 


6 


z 2 


2 


6, 6 


7 




3 


2, 2, 4 



TABLE IV: The irreps of Z 2 » D 6 . The number of an orbit 
corresponds to the number in Fig. |3] (.D&)^ is the Little 
Group of the orbit. 



D. Non-Abelian defect condensates in hexagonal 
crystals 



Let us now turn to the hexatic phase. Before turning 
to the Hopf symmetries involved in the next section, we 
first recall some of the standard analysis. The phase is 
two dimensional. It is translationally invariant, but the 
rotational group is broken to %§. Such a phase looks 
strange at first sight, because it seems that to break ro- 
tational symmetry we need to put the atoms on a lattice, 
thereby also breaking translational symmetry. We know 
that this is not the case because such a hexatic phase is 
actually realized in nature^. 

That it is not impossible to break translational sym- 
metry without breaking rotational symmetry, can be un- 
derstood from the representation theory^ of M 2 XI 50(2). 
A suitable order parameter field is provided by <fi — e l6e , 
where 9 measures the angle the bonds between neigh- 
boring molecules makes with the x-axis. In the high- 
temperature isotropic phase, < e t6e >= 0. In the ordered 
phase, < e t6e >^ because the angle of the bonds is on 
average a multiple of 2?. 

The hexatic phase can arise via spontaneous symmetry 
breaking from the isotropic phase, but it can also be the 
outcome of the melting of a two-dimensional hexagonal 
crystal. This is what happens in nature: as we increase 
the temperature of the crystal, at some point the trans- 
lational symmetry is restored, but locally the molecules 
are still arranged hexagonally. Of course, locally there 
is some translational symmetry, but the correlations of 
translations decay exponentially. 

We noted before that in a two-dimensional crystal the 
disclinations and the roton are very massive, and are not 
excited at low temperatures. As the temperature is in- 
creased and the translational symmetry is restored (sig- 
nalling the transition to the hexatic phase), the trans- 
lational Goldstone modes and the dislocations disappear 
from the spectrum. At the same time, the roton becomes 



massless, while the disclinations are the allowed defects 
in this phase. 

The traditional phase diagram of a two-dimensional crys- 
tal contains a line separating the crystal phase and the 
liquid phase. The phase transition is first order. How- 
ever, for the hexagonal crystal the transition to the liquid 
phase can occur via two second-order phase transitions: 
first a phase transition to the translationally invariant 
phase with broken rotational symmetry, and then a tran- 
sition to the liquid. 

The mechanism behind the phase transition from crys- 
tal to hexatic is well known, and quite remarkable: it is 
the formation of a dislocation condensate 14 . Defects in 
a phase with global symmetry breaking cost an infinite 
amount of energy, because their energy diverges logarith- 
mically with the size of the systemic. However, defect- 
antidefect pairs cost a finite amount of energy, and at 
finite temperatures they have a finite density, and an av- 
erage separation between the defects in a pair, which we 
call a. At the phase transition, a diverges, so that the 
single defects fill the ground state. This restores trans- 
lational symmetry. In the next section we will discuss 
why it restores symmetries. Here we merely note what 
happens next if the ground state is now filled with single 
dislocations. We've seen that these dislocations in prin- 
ciple can decay into a disclination-antidisclination pair. 
They also have an average separation, b. At the sec- 
ond phase transition, the disclinations blow out, thereby 
restoring rotational symmetry. 

The hexatic phase has been experimentally observed 
in numerous three dimensional phases, such as smectics, 
dense solutions of DNA, flux arrays of High-T c supercon- 
ductors, superfluids 4 , etc. All these phases have some 
kind of layered order. The hexatic smectic is also called 
the smectic B phase. The bond-orientational order is 
within the layers of the smectic. For references on hex- 
atic phases, we refer to the book by Nelson^. 



E. The hexatic phase 

Let us now apply our methods to this situation of the 
hexatic phase as an example of defect mediated melting. 
We start with the Hopf symmetry of the hexagonal crys- 
tal: 

A = F{I? x (Z 6 x Z)) <g> C((Z 2 x D e ) x Z), (21) 



and its representations which were given in Table ILLTI 

Now we consider a class sum defect condensate of the 
basic dislocation (1,0) belonging to the class C( a ^ r 6k 
with k = (see Figure 0J: 

\M = |(i,o)> + |(o, i)) + K-i,i)) 

+|(-1,0)) + |(0,-1)) + |(1,-1)). (22) 

It is not hard to determine what the intermediate T r and 
unconfined U symmetries are in this case (the derivation 
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is carried out in a related papers): 

%. = F(Z 2 x (Z 6 x Z)/Z 2 ) <8)C((Z 2 xi D 6 ) x Z) 
U = F{Z 6 xZ)»C(fl 6 xZ). (23) 

Zg x Z corresponds to disclinations, and Z?6 X Z corre- 
sponds to rotations and inversions. The Z part of Dq x Z 
corresponds to multiples of the 2n rotation. 

FIG. 4: An artist impression of the class sum dislocation de- 
fect condensate |0o) = 0)) + |(0, + 1)) + |(-1, 0)) + 
|(0, —1)} + |(1, —1)) with hexagonal symmetry. 

The particular condensate of dislocations is invariant 
under rotations, hence disclinations are not confined. Ba- 
sically the conclusion of our analysis is that the transla- 
tional modes of the lattice have completely disappeared. 
There are no more dislocations (because they are con- 
densed), and the translational phonons are confined. As 
we discussed earlier, this signals the restoration of trans- 
lational symmetries, so the unconfined and residual sym- 
metry algebras of the hexatic phase are actually 

% = F(Z 2 xi (Z 6 x Z)/Z 2 ) ® C((R 2 x D 6 ) x Z) 
U = F(Z 6 xZ)»C(l 2 xi (D 6 xZ)). (24) 

We have thus succeeded in interpreting the phase tran- 
sition from the hexagonal crystal to the hexatic phase (in- 
cluding also the translational symmetry restoration), as 
a Hopf symmetry breaking phenomenon through a par- 
ticular defect class condensate. 



under ISO{2) = R 2 xi SO{2). One can show that the 
two dimensional hexagonal and isotropic crystals have 
exactly the same Lagrangians 1 ? (they only have two in- 
dependent elastic constants), so we could consider this 
isotropic crystal to be the continuum limit of a hexago- 
nal crystal. 

Now the topological nematic is a phase obtained by a 
defect condensate in the isotropic crystal. To obtain a de- 
fect condensate that fully restores the rotational symme- 
try group 50(2), we must first realize that in the contin- 
uum limit the dislocations carry a continuous label. The 
Burgers vector of a dislocation in the continuum limit is 
characterized by a two-dimensional vector \(a,b)}, with 
a, b G R 2 . Thus a dislocation condensate that leads to 
restoration of the full rotational symmetry group is given 
by 

\4> r ) = dO \(cos6,sin6)). (25) 
Jo 

Strictly speaking, this condensate is the continuum limit 
of the class sum defect condensate we considered in the 
hexatic phase. It is indeed a class sum defect conden- 
sate of the isotropic crystal. It is in that sense that our 
analysis naturally incorporates the topological nematic. 

Furthermore, our analysis demonstrates how the par- 
tial restoration of the symmetry group can be understood 
by considering defect condensates that are not a class 
sum, as we will show next. Let us return to the hexatic 
phase and analyze the transition to the isotropic phase. 
The hexatic phase is a phase described by the Hopf alge- 
bra (JUJl 

A = F(Z 6 x Z) ® C(M 2 x (D 6 x Z)). 

Now assume a disclination condensate of |r), the 2l ro- 
tation, then the residual symmetries are successively 

% = F(e) <g> C(R 2 x (Z 6 x Z)) 

U = CR 2 . (26) 

G. A quantum smectic? 



F. A topological nematic 

We now wish to comment on a phase called the topo- 
logical nematic^. This is a two dimensional phase ob- 
tained from an isotropic crystal, where dislocations have 
condensed, but rotational symmetry is completely unbro- 
ken. In the defect condensate \\T2V that led to the hex- 
atic phase, the residual rotational symmetry group was 
Zg, but in the topological nematic the rotational sym- 
metry group is U{1). Clearly, an isotropic crystal only 
makes sense in the continuum limit, because by definition 
a crystal is not isotropic. However, in the continuum 
limit one can definitely have a Lagrangian of a crystal 
in terms of the displacement field u(x) , that is invariant 



We mentioned that with non-Abelian defects there are 
in general quite a few possibilities for the condensate. 
Our formalism allows us to determine what the physical 
implications are. We conclude by considering, instead of 
the class condensate (I22|l . the defect-antidefect conden- 
sate, 

|0„) = |(1,0)) + |(-1,0)}. 

The residual and unconfined symmetry algebras turn out 
to be 

% = F((Z 2 x Z 6 ) x Z)/0 x Z) <g> C(Z 2 x (Z 2 x Z 2 )) 
U = F((Z x Z 2 ) x Z) <g> C(Z x (Z 2 x Z 2 )). (27) 
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The dislocations along the (1, 0) direction are condensed, 
and the phonons with momentum in that direction are 
confined. Consequently the symmetry is restored in that 
direction, so that 

% = F((Z 2 x Z 6 ) x Z)/0 x Z) ® C(Z 2 x (Z 2 x Z 2 )) 
U = F((Z x Z 2 ) x Z) g) C((K x Z) x (Z 2 x Z 2 )). (28) 

This is a 2D smectic. Note that in a conventional clas- 
sical smectic, the Goldstone modes corresponding to os- 
cillations of the rods within a plane are massive (because 
of the analog of the Higgs phenomenon^). Determina- 
tion of whether the modes are massive or not requires a 
dynamical analysis. What we do know from our topo- 
logical analysis is that these Goldstone modes are not 
confined. 

The mass term for the Goldstone modes in classical 
smectics arises because the internal symmetry of the rods 
is coupled to the external symmetries (the Lagrangian is 
only invariant under coupled internal and external ro- 
tations). In the 2D smectic we have uncovered now, 
there are no internal symmetries. However, the defect- 
antidefect pairs in our condensate behave like rods in a 
conventional smectic. Note that they are coupled to ex- 
ternal symmetries, in fact they arise due to the break- 
ing of external symmetries. Thus we also expect the 
Goldstone modes associated with the oscillation of these 
"rods" to be massive, and only the longitudinal oscilla- 
tion of the planes to be massless. And so our analysis 
has naturally led us to the quantum smectic discussed 
by Zaanen and Nussino\»i&. 



IV. CONCLUSIONS 

In this paper we have applied the formalism of quan- 
tum double symmetry breaking to a number of physical 
phenomena linked to the formation defect condensates. 
We showed that both for cases with internal and exter- 
nal symmetries, nontrivial known and conjectured results 
can be recovered. It is especially fruitful to apply the 
method in cases where one is dealing with non-Abelian 
symmetries and consequently with non-Abelian defects, 
because in such situations many inequivalent conden- 
sates can be considered, which lead to effective theories 
with very different low energy degrees of freedom. It is 
also nice that the criterion for confinement determines to 
what extent a restoration of the original gauge symmetry 
takes place. The method enables one also to systemat- 
ically analyze dyonic condensates, but in this paper we 
have restricted ourselves to pure defect condensates. 

It is clear that there are many two dimensional (quan- 
tum) systems that can be systematically studied from 
this general point of view. In a paper appearing in paral- 
lel we treat non-Abelian nematic crystals with nontrivial 
point group symmetries of the tetrahedral, octahedral 
and icosahedral typei All systems in which conformal 
field theory plays a role, such as quantum Hall systems, or 



Bose-Einstein condensates, fall naturally in the category 
that should be analyzed. For example it would be in- 
teresting to investigate wether the plateau hierarchy pro- 
posed by HaldaneiS, where the plateaus are linked by the 
successive formation of quasi particle condensates, could 
possibly fit in our scheme. It may also be interesting to 
investigate the conceivable phases of (2+l)-dimensional 
gravity, along these lines. 
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Appendix: 
Inducing representations of 1? x De, 



Here, we illustrate the derivation of an irrep of Z 2 x 
Dq discussed in the main text, using Mackey's induction 
theorem (see any good book on representation theory^) . 
Take orbit nr 5 in fig. with momentum k = (■§,§) 
(i.e., (j-, ■£-) with a = 1). A rotation of 180 degrees 
(r ) gives an equivalent vector, i.e. a momentum that 
represents the same irrep of Z 2 . As a matter of fact, so 
do e, sr 2 , sr _1 , so £fri £ is isomorphic to Z 2 x Z 2 = (r 3 ) x 
(sr 2 ), where in general (g) is the group generated by g. It 
has four one-dimensional irreps. Now [Dq : Z 2 x Z 2 ] = 3, 
so we get four three-dimensional irreps. Denote an irrep 
of Z 2 x Z 2 by p m ,n, with m, n S Z 2 . Explicitly, a basis of 
the irreps is given by a set of representatives of left Z 2 x Z 2 
cosets, which we take to be {|e}, \r), \r 2 )} (the order is 
important when we write down the matrices explicitly). 
The matrices of the irrep are 



(a,b) 





V 



/ 
1 
V 1 



e 2 





(a+b) 




e 2 



(-a) 






bpif-(— a— b) 



e 2 "e 2 





f o 




S 1— > 









1 & mn 






As a representative example, let us illustrate the deriva- 
tion of the second column of the matrix corresponding to 
s (the derivations of the other columns is analogous): 



s\r) = \sr) = \r(sr 2 )) = \r)e l 



(29) 
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of 50(2). This leads to the conclusion that the irrep is 
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